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I. INTRODUCTION 



Understanding of the nonleptonic decays of light hadrons in the Standard Model (SM) 
has been one of the most difficult issues. While the SM explains the weak processes involving 
the change of strangeness by considering W exchange, the description of nonleptonic decays 
of light hadrons is rather complicated on account of the strong interactions in the low- 
energy regime (below 1 GeV). The problem is typified by the AT = 1/2 selection rule, 
best known as the fact that the isosinglet amplitude of the K ^ nn decay dominates over 
the T = 2 amplitude by about 22 times. Despite a great deal of effort this dominance 
of the AT = 1/2 channel over the AT = 3/2 it has not been explained in a satisfactory 
way. The effective weak Hamiltonian derived from evolving the simple ly-exchange-vertex 
from a scale of 80 GeV down to 1 GeV [§ 0, §, |, |, |, 0, |, | presents a part of the 
answer, showing that the perturbative gluons inherent in the Wilson coefficients enhance 
the AT =1/2 channel. However, the perturbative gluons alone are not enough to explain 
the AT = 1/2 rule completely. Thus, we need to consider other sources of the AT = 1/2 
enhancement. Since the structure of the light hadrons intimates already the importance 
of nonperturbative QCD, it is natural to consider its significance in describing low energy 
processes such as K ^ mi decay. 

Recently, we have investigated the effective AS* = 1,2 weak chiral Lagrangian to order 
0{p^) within the framework of the chiral quark model (xQM) UTOf, focusing on determining 
the low energy constants (LECs) in the effective weak chiral Lagrangian. However, the ratio 
of the gs/g27 obtained from the xQM deviated from the phenomeno logical values obtained 
in chiral perturbation theory (xPT). There was a suggestion to include the gluon condensate 
which is of order 0{asNc) |]11|. However, Ref. |]1T| did not consider the OlasNcY corrections 
which might be of significance at the same extent to the 0{asNc) corrections. 



In the present work, we want to improve our former study [|T0|, based on the more 



general effective low-energy QCD partition function derived from the instanton vacuum 
which pertains to nonperturbative QCD. The main feature lies in the fact that the coupling 
strength of the chiral symmetric quark- Goldstone interactions is momentum-dependent. In 
fact, switching off the momentum dependence of the coupling and adding an appropriate 
regularization scheme will lead to the usual xQM used in the former investigation . Due 
to the complexity of the formalism in the present approach we first concentrate on order 
C»(p2) and a part of 0{1/N^). 

As far as the large Nc expansion is concerned, while it explains the strong-interaction 
sector quantitatively well, the nonleptonic weak interactions defy any explanation from 
the strict large Nc limit. Because of the fact that the AT = 1/2 enhancement is badly 
underestimated in lowest oder in the large Nc expansion, it is inevitable to go beyond 
the leading order (LO) in l/Nc [|^, ^ ^ ^ ^ Moreover, if one expects a large 



contribution from the next-to- leading order (NLO) in 1/Nc, a convergence problem for the 



1/Nc expansion may occur p6[. In this case, one has to consider higher order corrections in 
l/Nc- Furthermore, since there are many different sources of nonleading-order corrections in 
1/Nc, one has to carefully analyze 0{1/Nc) corrections. However, dealing with higher order 
corrections and all contributions of 0{1/Nc) is beyond our work. Thus, we will restrict 
ourselves in this work to a part of the l/Nc corrections: The l/Nc diagrams from the quark 
operators. 

The instanton vacuum elucidates one of the most important low-energy properties of 



QCD, i.e. the mechanism of spontaneous breaking of chiral symmetry |T2|, O, O]. The 
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Banks-Casher relation tells that the spectral density z/(A) of the Dirac operator at zero 
modes is proportional to the chiral condensate known as an order parameter of spontaneous 
breaking of chiral symmetry: 

(^w ^ (1) 

The picture of the instanton vacuum provides a good realization of spontaneous breaking of 
chiral symmetry. A finite density of instantons and antiinstantons produces the nonvanishing 
value of z^(0), which triggers the mechanism of chiral symmetry breaking. The Euclidean 
quark propagator in the instanton vacuum acquires the following form with a momentum- 
dependent quark mass generated dynamically, identified with the coupling strength between 
quarks and Goldstone bosons: 

_ ^ + ^M{k) 

^^^^^ - WTmW) 

with 



M{k) = const- ^^^F\kp). (3) 

The ratio N/V denotes the instanton density at equilibrium and the p is the average size of 
the instanton. The form factor F{kp) is related to the Fourier transform of the would-be zero 
fermion mode of individual instantons. As will be shown later, the momentum dependence of 



the constituent quark plays an essential role in improving our previous results [10| concerning 



the effective weak chiral Lagrangian. It is in line with the recent works on the pion wave 



functions [T^, |T6| and skewed parton distribution [|Tj], where the momentum-dependent 
quark mass plays a crucial role as well. 



The instanton vacuum induces effective 2A^j-fermion interactions [0, 0. For example. 



it has a type of the Nambu-Jona-Lasinio model for Nj = 2 while for Nj = 3 it exhibits 



the 't Hooft determinant [jT9|. Goldstone bosons appear as collective excitations by quark 
loops generating a dynamic quark mass. Eventually it is found that at low energies QCD 
is reduced to an interacting quark- Goldstone boson theory given by the following Euclidean 
partition function |^ 

' ^ ' (4) 



where If^^ stands for the pseudo-Goldstone boson: 

U^^ix) = Uix)^^ + U^(x)^^ = exp (m'^(x)A'^75) • (5) 

The a is the color index, a = 1, - ■ ■ ,Nc and / and g are flavor indices. M{—id) is the 
constituent quark mass being now momentum-dependent. It will play a main role in the 
present work. This effective theory of quarks and light Goldstone mesons applies to quark 
momenta up to the inverse size of the instanton, p~^ ~ 600 MeV, which may act as a scale 
of the model (PxQm)- A merit to derive the xQM from the instanton vacuum lies in the 
fact that the scale of the model is naturally determined by p~^. Furthermore, mesons and 
baryons can be treated on the same footing in the xQM- Fo'^ example, the model has been 
very successful in describing the properties of the baryons ||21||. 
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We want to mention that there is a problem related to the gauge invariance in Eq.(^). 
Due to the nonlocality of the interaction expressed in Eq.(Q), the vector and axial- vector 
Ward-Takahashi identities are not satisfied, that is, the conservation of the vector current 
(CVC) and the partial conservation of the axial- vector current (PCAC) do not hold. Thus, 
one has to make the effective action gauge-invariant. A well-known prescription to solve this 
problem is to insert the path-ordered exponent into the effective action [53, 123, RDI, pi], B^l, 



At the scale of the W boson, the charged- current weak interaction of the quarks is mediated 
by the W boson only and thus is a product of two local currents. However, when the 
charged- current weak interaction is scaled down to the hadronic scale, i.e. /i = 1 GeV, the 
renormalization of the QCD quantum corrections make the current-current operator mixed 
with other different types of twist-four local composite operators However, in the strict 
large limit, the contribution of all other composite operators is suppressed by l/N^ except 
for the current-current one which is the same as that at the W scale. Thus, it need not be 
renormalized and the conservation of the currents must be taken into account. 

The outline of the present paper is as follows: In Section 2 we sketch the instanton- 
induced chiral quark model, emphasizing in particular the momentum dependence of the 
constituent quark mass and explain how to perform the derivative expansion in the presence 
of the momentum- dependent constituent quark mass. In section 3, we show the basic for- 
malism to obtain the effective weak chiral Lagrangian starting from the effective weak quark 
Hamiltonian. In section 4, we dicuss the large Nc limit in the nonleptonic weak Hamilto- 
nian and the conserved currents with the nonlocal interaction. Section 5 is devoted to the 
derivation of the effective weak chiral Lagrangian to order 0{p'^) and (9(l/iVc). In Section 
6 we discuss results. The conclusions and outlook are given in Section 7. 



II. EFFECTIVE CHIRAL THEORY FROM THE INSTANTON VACUUM 

The effective low-energy QCD partition function in Euclidean space can be written as 



The M{k) is the momentum-dependent constituent quark mass expressed as follows: 

M{k) = MoF\k/A). (7) 

If we choose the F{k/A) to be constant and add a regularization {e.g. Pauli-Villars or proper- 
time), the partition function becomes just that of the usual xQM. The original expression 
for the F{k/A) ||T2| , p!3| , which is obtained from the Fourier transformation of the would-be 



zero fermion mode of individual instantons with the sharp instanton distribution assumed, 
is as follows: 

F{k/A) = 2z [lo{z)K,{z) - h{z)Ko{z) - -J,{z)K,{z)^ . (8) 

Here Jq, /i, Kq, and Ki denote the modified Bessel functions, z is defined as 2; = k/(2A) 
and the cutoff parameter A is in this case just the inverse of p. When k goes to infinity, 
F{k/A) = F{kp) becomes 
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Actually, the momentum-dependent quark mass is related to the nonlocal regularization 
in Euclidean space. There are other ways of understanding the nonlocal effective interaction 
without relying on the instanton vacuum ||3^ . In those methods, the momentum-dependent 
quark mass as a regularization appears as delocalizing the quark fields. So, various types 
of the M{k) as a regulator with the regularization parameter A has been used by different 
authors [^, |T^, . 

Therefore, we will not confine ourselves to the expression given in Eq.(|^) but rather try 
three different types of the M{k): 



M{k) 



Eqs.(7,|) 

M A* 

Moexp (-^ 



(10) 



The M{k) is normalized to Mq aX k = 0. Originally, M(0) is taken to be around 350 
MeV corresponding to the values of the following parameters: ^ ^ 1 fm and p ^ 0.35 fm. 
The R is the average distance between neighboring instantons. However, we will regard Mq 
as a free parameter ranging from 200 MeV to 450 MeV and fit for each Mq the parameter 
A to the pion decay constant = 93 MeV. Figure 1 shows the momentum dependence of 
the three different types of M{k) with Mq = 350 MeV. 

Having integrated out the quark fields in the partition function, we obtain 



Z = J P7r'^exp(-5effM), 



where the S'efrfvr''] stands for the effective chiral action: 

S-effH = -NJndetDiU'''). 
Here, the D{U"'^) is the Dirac operator defined by 

D = ip + i^M{-id)W'^M{-id). 



(11) 



(12) 



(13) 



The Dirac operator is not Hermitian, so that it is useful to divide the effective action into 
the real and imaginary parts: 



ReS'eflf 
ilmS'eff 



(S'eflf + SIq) 



efT 



eflfJ 



--Ncln det 
2 

--A'clndet 
2 



D^D 
D/D 



(14) 



In order to calculate the real part, we substract the vacuum part and use the derivative 
expansion. We therefore write 



Re^effkl -Re5efr[0] 



.^Trlnf^^l 
r f d^k 



2 

' 2 



d'^x 



ikx 



(27r)4 
d^k 



tr In 



Akx 



\Dl{d^d + ik)Do{d^d + ik)J 
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Here we used a complete set of plane waves for the calculation of the functional trace, 
summed over all states and took the trace in x. 'tr' then denotes the usual matrix trace 
over flavor and Dirac indices. The rhs. of Eq.(p!5[) can now be expanded in powers of the 
derivatives of the pseudo-Goldstone boson fields, , and of 2ik ■ d + d"^. After some 
manipulation using 



D\d ^ d + ik)D{d ^ d + ik) = (ip - ft - iJM{-id + k)U-^'JM{-id + k) 



x(ip-lt + i^M{-id + k)U^^ ^M{-id + k) 



-d' + e -2ik-d + JM{-id + k)U-^'M{-td + k)U^'^JM{-id + k) 



M{-id + k) ipU^') ^M{-id + k) 

d'^ + k'^-2ik-d + M^ -M {^W') - 2iMM'ka hda + U'^'" {dJJ^'' 



- MM' 



2d' + U-^" [d'W' + 2U-^^ (daU^') da 



2MM"kak^ 
2M'^k„k 



+ iM'k^ 



dI {d^d + ik)Do{d ^d + ik) = -9^ + k"^ + - 2ik ■ d 



-AiMM'k ■ d - 4:M'\akfsdadp - 2MM'd'^ - AMM"kakpdadp + 0{d'^ 



where 



M = M{kl M' = i^, M" 



1 



'd^M{k)^_dM{k)\^ 



dk^ 



dk 



(16) 



(17) 



(18) 



2k dk 

we obtain eventually the effective chiral action with the momentum-dependent quark mass: 

J n 



ReS-efffTT^ - ReS'efr[0] = J d^x {L^L^) . 



(19) 



In Eq.([T9|) () denotes the flavor trace, = iU'^dfjU is a Hermitian Nj x Nf matrix and /tt 
denotes the pion decay constant: 



J n 



d^k M^{k) - \M{k)M\k)k + lM'\k)e 



(2vr) 



with 



M'{k) 



{k^ + M^{k)y 

dM{k) 
dk ' 



(20) 



(21) 



When we switch off the momentum dependence of the constituent quark mass, we end up 
with the former expression for |[lO[| : 



J TT 



AN' 



d'^k 



(27r)4 (P + M'^y 



M 



const. 



(22) 



which is logarithmically divergent. The quark condensate, which is just the trace of the 
quark propagator given by Eq.(^, is written by 



d^k 



M{k) 



(27r)4F + M(/t)2' 



(23) 
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The subscripts M and E stand for Minkowski and Euclidean space, respectively. The gluon 



condensate is expressed as follows |13|: 



TT 



32iV, 



(27r)4p + M2(A;)' 



(24) 



It is already known that the imaginary part of the effective chiral action is identical to 
the Wess-Zumino-Witten action |3^ with the correct coefficient, which arises from the 
derivative expansion of the imaginary part to order 0{p^) An appreciable 

merit of using the momentum-dependent quark mass as a regulator was already pointed 



out by Ball and Ripka The momentum-dependent quark mass provides a consistent 
regularization of the effective action given in Eq.(|T^ in which its real and imaginary parts 
are treated on the same footing and thus pertinent observables such as anomalous decays 
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27 are safely recovered even if M{k) acts as a regulator. 



III. EFFECTIVE WEAK CHIRAL ACTION 

In this section, we will show how the effective weak Hamiltonian is incorporated into the 
present framework. The effective weak chiral action S'^'^~^'^[7r"] can be written as follows: 



exp = j Vi)Vi)^ exp j d^x {^iP^DiP - U, 



/A5=l,2 
'•cff 



(25) 



Here the effective weak quark Hamiltonian Tifs consists of ten four-quark operators : 

= -^K.K;Ec.(/^)Q.(/i)+h.c.. (26) 

The Gp is the well-known Fermi constant and Vij denote the Cabibbo-Kobayashi- 
Maskawa(CKM) matrix elements. The r is their ratio given by r = —VtdV^l/VudV*^. The 
Wilson coefficients Cj(/i) are defined as Cj(/i) = Zi(/i) +Tyi{ii). The functions Zi{fi) and yi{fj.) 
are the scale-dependent Wilson coefficients given at the scale of the /i. The Zj(/i) represent 
the CP-conserving part, while yi{fi) stand for the CP-violating one. The four-quark op- 
erators Qi contain the dynamic information of the weak transitions, being constructed by 
integrating out the vector bosons and Z and heavy quarks t, b, and c. The four-quark 
operators are given by 

Qi = -4 [sI^^Plu^) {u^p^^Pidc) , Q2 = -4 (sl^^PLUa) {ul^^Pidp) , (27) 

q=u,d,s q=u,d,s 

Q5 = -4 [si^^PLda] J2 [(lh^^Pmp) . Qe = -4 (si^i^PLdp) ^ {qll^.PRqo 

q=u,d,s q=u,d,s 

Q7 = -6 {sl^^^PLda) E {qlQl^^PRQp) , Qs = -6 (sh.PLdp) E {qIQi^PrQc 

q=u,d,s q=u,d,s 

Qg = -6 (sl^^PLda) E {QlQl^^PLqf3) , Qio = -6 (sb^PLdp) E {^iQlf^PLqc 

q=u,d,s q=u,d,s 

where Pl,r = | (1 i 75) are the chiral projection operators and Q = idiag(2, —1, —1) denote 
the quark charge matrix. The Qi and Q2 come from the current-current diagrams, while 
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Qa to Qe 0! ilj si^nd Q7 to Qio [§] are induced by QCD penguin and electroweak penguin 
diagrams, respectively. Note that only seven operators in Eq.(p7D are independent. For 
example, we can express Q4, Qg, and Qio as follows: 



Q4 



■Q1 + Q2 + Q3, Qg 



2 (3Qi - Qs) 



Q 
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Q2 + -(Qi 



Q3 



(28) 



Under the chiral transformation SU(3)l x SU(3)^ the four-quark operators Q2 — Qi, Q3,4,5,6 
transform like (8^, Ir). The Qi + Q2, Qg.io transform like the combination of (8^,, Iji) and 
(27l, 1r), while the Qj^s transform like {Sl,Sr). 



The AS* = 2 effective weak Hamiltonian is expressed as |45 



o/A5=2 



with 



= \lr]iS 



167r2 



7712 



T (Ac, At, m\, m\, M^) QA5=2(/i) + h.c. 



mi 



mr 



(29) 



(30) 



and the parameters Ag = Vg^V^t, denote the pertinent relations of the CKM matrix elements 
with q = u,c, t. The functions Si are the Inami-Lim functions |^, being obtained by 

integrating over electroweak loops and describing the \^S\ = 2 transition amplitude in the 
absence of strong interactions. The b{n) is again the corresponding scale-dependent Wilson 
coefficient. The coefficients rji represent the short- distance QCD corrections split off from 
the 6(/i) The four-quark operator Qas=2 is written as 



Q 



AS=2 



-4 [si-f^PLd^) {sli^PLd 



(31) 



Since the Fermi constant Gp is very small, one can expand Eq.(^) in powers of the Gp 
and keep the lowest order only. Then we can obtain the effective weak chiral Lagrangian 



AS=1,2 
eff 



1 

Z 



J 



(32) 



If you write a generic operator for the four-quark operator Qi for a given i in Euclidean 
space such as 

Qi{x) = -^^(x)7^P^j,iAi^(x)^t(a;)7^PH,LA2^(x), (33) 

where Ai 2 denote the flavor spin operators, then we can calculate the vacuum expectation 
value (VEV) of Qi{x) as follows: 



[Qi) = ^1 Vm^Qi{x)exp U d^z^lj^DiP 



1 



X exp 
r(i) , 



L 



d^k 

(2) 



Jk(x-y)_ 



64'\x)64'\y) 



ti\nDMz),J2{z)) 



Jl = J2=0 



Here, D is 



DiJi{z),J2{z)) = + J«(z)7„Pk,lAi + (z)7^Pk,lA2 



(34) 
(35) 
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The hf^ and hf^ are given by 



(1) 



(2) 



-Nhi 



N^tT 



X ) { X 



X 



0{K 



1,4,6,8,10 



X 

(^1 



D 



tr 



X 



D 



1^,Pr,l^2 



X 



+ O (iVj 



tr 



2,3,5,7,9 



(36) 



(37) 



where Ai^2 are the corresponding flavor matrices. Applying the derivative expansion to the 
operators (^1,2)^ up to order (9(9^) they can be written as 



{A 



1,2; 



^ - r / 1 \ 

X \2i (1 + 2MM') k^da + M{^W^) + 2iMM'kJJ-^'{dJJ^'- 
'\ + 2MM') 92 + MM' (f/"^5(92f/^5) + 2f/~^5((9„f/^5)(9„ 



+ 4(MM" + M<^)kjzfi {da_d(i + U-^^^{d^U^')dp + ^U-^' {d^d^U^^^ 



X 



+ i\M"-]^^\kMdAU-~<^' 



7m^l,rAi,2, 



(38) 



where M, M' and M" are the fc-dependent functions defined in Eq.(^). Before we proceed 
to determine the low energy constants for the effective weak chiral Lagrangian, we want to 
discuss the role of the large N^. limit and related problems of the current conservation. 



IV. THE l/Nc EXPANSION AND CURRENT CONSERVATION 



The large expansion enters into two different places: Wilson coefficients and four- 
quark operators. The dependence of the 6x6 anomalous dimension matrix 7 in the 
basis of Qi, 2,3,4,5,6 was shown by Ref. |25|. The operators Q7,8,9,io are not mixed with 
Qi, 2,3,4,5,6 05 since the LO and NLO anomalous dimensions 7ij, i = 1, ■ ■ ■ 6, j = 7, ■ ■ ■ 10 are 
zero. Up to (9(1/A^c); we write the matrix of the anomalous dimension [| 



51 2e 



as 



with 



/O 3 
3 

3 


\^0 





1 

A 



3 


3 



1 tts 

2 TT 










(39) 





1 



3 

-3 



ilL (^-3A^, 



(40) 
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where as(~ 7^) denotes the strong running couphng constant. In the strict large- A''^ hmit, 
all elements vanish except the anomalous dimension of the Qq operator 756. It means that 
the operator Q2 is not mixed with the other operators Qi{i 7^ 2) and Qq is only renormalized. 
Because of ~ The Wilson coefficients Ciii 7^ 2) go to zero and C2 becomes unity. Thus, 
in the strict large-A^^c limit the operator Q2 does not require any renormalization, so that 
it restores the original current-current form of quark charged weak interaction at /i = Mw, 
i. e. the factorized form. 

Strictly speaking, the Q2 is the only contribution to the effective weak chiral Lagrangian 
in the large A^^^ limit and other operators must be treated as the NLO corrections according to 
the Wilson coefficients ||26| . However, we will take the weak Hamiltonian as our starting point 



more practically without considering the large Nc behavior of the anomalous dimensions 

Since the momentum- dependent quark mass presents the nonlocal interaction between 
the Goldstone-bosons and quarks, the vector and axial-vector currents are known to be 
not conserved. Making the effective action gauge-invariant with some approximations, one 



is able to derive conserved currents [30, pl^, 32, B3l. The conserved vector and axial 



vector currents in Euclidean space with the momentum-dependent quark mass are written 
as follows p3| : 



where yM^ denotes dyM/dp^. Last two terms in and are required so that the 
currents can be conserved. 

The pion decay constant /,r is related to the following transition matrix elements: 

{Q\AI{x)\tx\p)) = ^^J9^e^P■^^'^^ (42) 

where ^^s defined in Eq. f^ll) . The obtained from Eq.(|4^) is exactly the same as 
Eq.(pOl) p3| , which indicates that the PCAC is well satisfied with the additional nonlocal 



term in Eq. (]41|) . With the nonlocal terms turned off, we would end up with the Pagels- 
Stokar condition for fl 

' J (27r)4 (p + M2)2 ^ ^ 



Although the pion decay constant given in Eq. (|20|) is the correct one with the momentum- 
dependent quark mass, we want to use the Pagels-Stokar condition for the normalization of 
the effective chiral Lagrangian for convenience. The reason lies in the fact that by using the 
Pagels-Stokar condition we need not consider the additional nonlocal part of the currents 
in deriving the VEV of the quark operator, since we obtain the same results as we use the 



conserved currents given in Eq.(pD and hence we reproduce precisely the large Nc results 



for the Bk factor as well as the LECs. 

V. EFFECTIVE WEAK CHIRAL LAGRANGIAN 
A. Leading order in the l/N^ expansion 

In order of 0{p^), only the operator Qg has the nonvanishing term: 

(Q8 + Q^)o(p0) = AmlM\U\^U^Q). (44) 
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However, since we are mainly interested in the LECs gg and 5^27 of the effective weak chiral 
Lagrangian in LO, we proceed to calculate the NLO terms, i.e. those in 0{p'^) order: 



+ k;i/ci(p2) 



{\,L,L,) + ^f^{\lL,){XiL,) 



0(p2) 

+ y7/ci(p2) 

0(p2) 



{Qs + Ql) 



\yio + yio/c)(p2) 



AS'=2 



0{Nl) 
A5=2/Ci(p2) 



0, 

16iV2/C2(A6L^L^), 
0, 

UNlM {V + n) (AeL^L^), 

24iV2/C2(L^A6)(i?^g), 

ASNlM {V + 7^) [(f/Ae {d^U^) (d^U) U^Q) 

IQNjie {^-{KL,L,) + lt}t(A|^,)(ALL,) 
16iV2/c2(A6L^A6L^), 



(45) 

(46) 

(47) 
(48) 
(49) 
(50) 
(51) 

(52) 
(53) 

(54) 
(55) 



where f-i are the eikosiheptaplet projection operators |10, 4£] and the functions /C, A^, P, 
and TZ are expressed as follows: 



7^ 



(27r)4 (A;2 + M2)' 
(27r)4 F + M2 ^ 



/.'(PS), 



4A^. 



oZlvl I 4 I z 2 4 



(2^ 



P + M2 (p + M^y 



{k^ + M^y 



(56) 



(27r)4 (A;2 + M2)^ 

Note that the function /C is just the same as the Pagels-Stokar condition and the function 7W 



is directly related to the quark condensate. The AS" = 1 effective weak chiral Lagrangian ||50 
is given in Minkowski space as follows: 



^A5=l,C'(p2) 



-eff 



Gf^. .4 / (i/2)^n/2i 



(1/2) 
8 



(1/2)^(1/2) (3/2)^(3/2) 
y27 ''-27 "T i/27 ''-27 



(57) 



where 



XlL^Lf") + h.c. 
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iTfl 



h.c. 



-C 



(1/2) 
27 



5 ^(3/2) 



(58) 



and 



^(1/2) 

A3/2) 
*-27 



XlL^L") + h.c. 

2 



A?L,)(a^L'^) + 2(a^L,)(a}L'^ 



5 (^IL^ 



A^L^\ + h.c., 



A^L'^) +h.c. 



(59) 



The Qs and (727 are dimensionless LECs of which the numerical values can be extracted from 
the CP-conserving K —>■ mc decay rate and the AT = 1/2 enhancement is reflected in these 



constants. From the analysis in chiral perturbation theory with chiral loops considered |51 



52| , the LECs have the following values 

IS'slcxp ^ 3.6, \g27 



cxp 



0.29, 



1^8 1 



cxp 



lfl'27 



12.5, 



(60) 



exp 



Comparing Eq.(|^ with Eqs.( P5| - plD , we determine the values of the LECs gs and g27 to 
the LO in the 1/Nc expansion: 



2 3 3 2 QAN^MiP + n) 

98 = -^Ci + -C2 + C4 - -Cg + -Clo H Jl Cq, 

3 3 9 9 

Q27 = -C^ H Co H Co H Cin. 

- 5 5 10 10 



(61) 



In the strict large Nc limit in which the Wilson coefficient C2 survives only and becomes one, 



we correctly reproduce the following large N^. results 1^1 

3 



fi'27 



Nc^oo 



Nc^oo 



5' 
3 



The effective AS* = 2 weak chiral Lagrangian to order 0{p^) is derived as: 



AS=2,0(p^) 
eff 



w 



47r2 



^ (Ac, A 



(62) 



(63) 



where Bk is known as the scale-independent Bk factor which is defined as Bk = BK{fJ')b{fi). 
The scale-dependent Bk{^ji) is related to the matrix element for the — mixing 



{K^\Qas=2\K^) 
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BK{fJ')fKmK, 



(64) 



which governs the — mixing at short distances. Here, fx and mx are the mass and 
decay constant of the neutral kaon, respectively. Comparing Eq.(^) with Eq.( |55|) in the 
chiral limit (/^r = /i<-), we immediately obtain the B^ factor in the large N^. limit |54 



B 



K 



3 
4' 



(65) 



12 



B. 0{Nc) corrections 

Now, we proceed to add the next-to-order corrections in the large Nc expansion. Though 
there are many different origins of the NLO corrections in the Nc expansion, we are not 
able to consider all possible NLO corrections. Here, we will restrict ourselves the 1/Nc 
corrections from the quark-quark operators. The VEVs of the quark operators in the NLO 
are as follows: 



1^4 + y4/o(p2) 
\t£^6 + t^6/0(p2) 



^^9 + ^9/0(^2) 



1 



= IQNcK' [^{X,W + -tf,{XlL,){XiL,) 

= 0, 

= Q4NcM{V + n){XeL^L^), 

= 0, 

= 48NcM {V + 7^) [(f/Ae (d^U^) (d^U) U^Q) 

+ {{d^U)(d,U^)uX,U^Q)], 



(Qio + Qlo) 



AS=2 



0{N^) 
A5=2/0(p2) 



+ Qa5=2) 



16iV,/C2 ^--{\,L,L,) + -t;.t(AiL,)(AiL,) 
16iV,/C2(A6L^A6L^). 



(66) 

(67) 

(68) 
(69) 
(70) 
(71) 

(72) 
(73) 
(74) 

(75) 
(76) 



The LECs and g27 are then obtained as follows: 



2 3 \ /3 
H — C5 4 



5iV, 



C2 + TT^S + C4 + 



J TT 



J TT 



C65 



527 



1 



1 + TT 7C1 + 7C2 + -7:^9 + tttCio 



10 



10 



3 2 

"5 ^ 5iVc 



C9+ 7 



5iV, 



ClO 



(77) 



If we turn off the momentum dependence of the M{k) and introduce a regularization with 
the cut-off parameter A for the quark-loop integrals, we end up with the former results in 
Ref. m-. 



/C 
M 
V 



J TT 

Wc 

J TT 



+ 



M 



8NrM 647r2 
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7^ 



M 
'32^ 



(78) 



Thus, the LECs ( ff?]) are reduced to those with the constant constituent quark mass |T0 



_2 

~5 ^JNc 



3 



+ 



5 SA^c 
C5 -|- 



C2 + J^Cs + C4 + 



4/4^2 



3 2 

C6, 



2 3 



5^27 



1 + 



3 3 9 9 

-Ci H Co H Co H Cin I . 

5 5 10 10 ' 



ClO 



(79) 



From Eq.(^), the factor with the 0{Nc) correction becomes 



B 



K 



(80) 



which is just one with Nc = 3 |E5l and agrees with the result of l/Nc approach [El, ISH . 



VI. RESULTS AND DISCUSSION 

We employed the Wilson coefficients Cj obtained by Buchalla et al. ^ . There are three dif- 
ferent renormalization schemes in Ref. For our best fit, we choose the naive dimensional 
regularization (NDR) scheme in this work. 

While in the former calculation based on the usual chiral quark model 0] three pa- 



rameters are involved, namely the pion decay constant, the quark condensate, and the 
constituent quark mass, the present results given by Eqs. (|6T| , [77|) include four functions /C, 
A^, P, and TZ. However, those functions depend on the M{k) containing the A and Mq, 
so that we have only two free parameters. Furthermore, the cut-off parameter A is fixed 
by reproducing the pion decay constant in our calculation. The value of Mq is more or less 
constrained by demanding the calculated quark condensate ( p3D and gluon condensate 
to lie inside the empirical hmits of -(350 MeV)^ < (ipi!) < -(200 MeV)^ [0, |8|, |5|] and 



(350 MeV)^ < i^GG) < (400 MeV)^ |60|, g, |2[. Thus, we first want to examine the 
dependence of the quark and gluon condensates on the momentum-dependent quark mass 
(see Eqs.(^,gl). 

In Figure 2 we show the dependence of the quark condensate on the Mq with several 
different types of the M{k) given in Eq.(|T^. The M{k) by Diakonov and Petrov [^, [TB 



and that of the dipole type produce the very similar results of the {il^i/j) pBI, the Gaussian 



type brings out noticeably smaller value of the quark condensate than the other two. From 
these results one can easily see that the M[k) with stronger tail produces the larger value of 
the quark condensate. The shaded band depicts the empirical range of the quark condensate. 

In Figure 3 we draw the dependence of the gluon condensate on the Mq in the same 
manner. As in the case of the quark condensate, the dependence on the different types of 
the M{k) is similar. Again, the shaded band represents the empirical range of the gluon 
condensate. Thus, examining the dependence of the quark and gluon condensates, one 
constrains the range of the value of Mq. 

In Figures 4 and 5 we present dependence of the LECs gg and (727 on the dynamical quark 

mass, with the NDR scheme employed {A.^^^ = 435 MeV). As seen in Eqs. (|6T| , [77D , the (727 
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does not depend on Mq. On the other hand, the gs shows a strong dependence on the Mq. 
It is due to the fact that the penguin operator contributes only to the gs and it brings about 
three functions containing the Mq and the quark condensate. The behavior of the gs with 
different types of M{k) is similar to the quark and gluon condensates. In particular, the gs 
goes up drastically below Mq = 250 MeV. As a result, the ratio gs/g27 exhibits a similar 
dependence on the Mq as shown in Fig. 6. 

Figures 7 and 8 depict the contribution of the NLO corrections 0{Nc) to the LECs gs and 
g27. The 0{Nc) contribution to the gs is very tiny, while it does to the (727 almost 30% and 
thus it makes the g2T_ deviate from the empirical value. As already discussed in Ref. |TD[ , the 
0{Nc) corrections make the ratio gs/g27 worse as shown in Fig. 6, compared to the empirical 
value (dotted line). We only can reach the empirical value below around Mq = 230 MeV 
which is not preferable according to the proper range of the Mq from the quark and gluon 
condensates. 

Now, we want to analyze the contributions to the LECs from the different quark operators 
Qi given in Eq. p7|) . These considerations are basically independent of the actual form of 
M{k). Since the penguin operator Qg contributes only to the octet part of the coupling, 
it contributes directly to the AT = 1/2 enhancement. In fact, the appearance of this 
octet penguin operator indicates already that the effect of perturbative gluons is of utmost 
importance in explaining the AT = 1/2 rule. Figure 9 shows each contribution of the quark 
operators Qi, starting from the Q2 which has the same form of the bare operator at the 
scale of the M^y. The operator Q2 has the largest Wilson coefficient in any renormalization 
scheme. The Qi takes the second largest Wilson coefficient. As drawn in Fig. 9, the 
LEG gs is much underestimated with the operators Qi and Q2 only and furthermore they 
are stable to the change of the Mq. Adding the penguin operator Qq brings the gs up 
dramatically and its dependence on the Mq is also very noticeable. The other contributions 
are negligibly tiny. Hence, the penguin operator Qg plays an essential role in enhancing the 
octet coupling. Actually, considering the fact that the M{k) is pertinent to the zero mode of 
individual instantons, the strong dependence of the penguin operator on the M indicates the 
importance of the effect of a part of nonperturbative gluons in describing the nonleptonic 
kaon decays. 

In the case of the (727, the operators Qi and Q2 are dominant and the Wilson coefficients 
Cg and Cio are negligibly small. 



VII. SUMMARY AND CONCLUSIONS 

In the present work we have investigated the effective weak chiral Lagrangian for AS* = 
1, 2 concentrating on the calculation of gs and (?27- We used the chiral quark model from 
the instanton vacuum as framework and incorporated the weak interaction by the effective 
Hamiltonian of Buchalla, Buras, and Lautenbacher The calculation has been done in 
a first step to order 0{p'^) and to LO and NLO in the l/N^ expansion. In contrast to the 
previous work by Ref. [T^, we used a momentum-dependent constituent quark mass M{k) as 
it arises from the instanton vacuum, which makes gs be closer to its corresponding empirical 
value. However, the g27 is untouched by the present work, so that it is not at all improved, 
since it is independent of M. 

For their ratio gs/ 927, this simple feature enlarged the values by about 50 % to 100 % and 
hence very much improved the theoretical values compared with the empirical ones. In fact 
we tried various shapes of M{k). However, the best values were obtained for the M{k) from 



15 



the instanton model of Diakonov and Petrov [|T2|, ^ and Mq ~ 230 MeV. Altogether the 
empirical value of gs/ g-ri with the Mq constrained to the quark and gluon condensates is still 
underestimated and one is still away from an explanation of the AT = 1/2 rule. However, 
the conclusion is very clear: If one wishes to derive the effective weak Lagrangian from the 
instanton vacuum of QCD by means of the chiral quark model the use of a momentum- 
dependent quark mass seems to be indispensable. 
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Fig.l: Momentum dependence of the constituent quark mass M{k). The sohd curve denotes 
the original M{k) — MoF'^{k/A) by Diakonov and Petrov, while the long-dashed one stands 
for the dipole-type M (k) and the dot-dashed one depicts the Gaussian one. The A has been 
fitted to reproduce = 93 MeV. 
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Fig.2: Quark condensate as a function of M[k). The solid curve denotes the original 
M{k) — MoF'^{k/ A) by Diakonov and Petrov, while the long-dashed one stands for the 

dipole-type M{k) and the dot-dashed one depicts the Gaussian one. The shaded band 
designates the range of the empirical values of the quark condensate. 
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Fig.3: Gluon condensate as a function of M{k). The solid curve denotes the original 
M{k) — MoF'^{k/ A) by Diakonov and Petrov, while the long-dashed one stands for the 

dipole-type M{k) and the dot-dashed one depicts the Gaussian one. The shaded band 
designates the range of the empirical values of the gluon condensate. 
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Fig.4: The low energy constant gs as a function of M{k). The sohd curve denotes the 
original M{k) — MoF^(fc/A) by Diakonov and Petrov, while the long-dashed one stands for 
the dipolc-typc M{k) and the dot-dashed one is the Gaussian one. The empirical value of 
I gig I is approximately 3.6 which is drawn in the dotted line. 
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Fig.5: The low energy constant as a function of M[k). It is independent of the M[k). 
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Fig.6: The ratio g%/ g-ri as a function of M{k). The sohd curve denotes the original 
M{k) — MoF^(/c/A) by Diakonov and Petrov, while the long-dashed one stands for the 
dipole-type M{k) and the dot-dashed one is the Gaussian one. The empirical value of 
Is's/S'^T.I is approximately 12.5 which is drawn in the dotted line. 
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Fig. 7 The LO and NLO contributions to the gg in the large Nc expansion. The long-dashed 
curve represents the LO contribution, the dot-dashed one the NLO contribution, and the 
solid one the full result. The empirical value of Igg] is approximately 3.6 which is drawn in 
the dotted line. 
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Fig. 8 The LO and NLO contributions to the g-ri in the large Nc expansion. The long- 
dashed curve represents the LO contribution, the dot-dashed one the NLO contribution, 
and the solid one the full result. The empirical value of \g^ is approximately 0.29 which is 
drawn in the dotted line. 
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Fig. 9: Contributions of each quark operator Qi to the low energy constant g^g as a function 
of M{k). The long-dashed curve denotes the contribution of the Q2, while the short-dashed 

one draws that of the Qi + Q2- The dot-dashed one depicts the with added and the 
solid curve represents the full result. The empirical value of l^fgl is approximately 3.6 which 
is drawn in the dotted line. 
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